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ABSTRACT 

Observations of hot Jupiters around solar-type stars with very short orbital periods 
(~ 1 day) suggest that tidal dissipation in such stars is not too efficient so that 
these planets can survive against rapid orbital decay. This is consistent with recent 
theoretical works, which indicate that the tidal quality factor, Q^, of planet- hosting 
stars can indeed be much larger than the values inferred from the circularization 
of stellar binaries. On the other hand, recent measurements of Rossiter-McLaughlin 
effects in transiting hot Jupiter systems not only reveal that many such systems have 
misaligned stellar spin with respect to the orbital angular momentum axis, but also 
show that systems with cooler host stars tend to have aligned spin and orbital axes. 
Winn et al. suggested that this obliquity - temperature correlation may be explained 
by efficient damping of stellar obliquity due to tidal dissipation in the convection zone 
of the star. This explanation, however, is in apparent contradiction with the survival of 
these short-period hot Jupiters. We show that in the solar-type parent stars of close-in 
exoplanetary systems, the effective tidal Q^, governing the damping of stellar obliquity 
can be much smaller than that governing orbital decay. This is because for misaligned 
systems, the tidal potential contains a Fourier component with frequency equal to the 
stellar spin frequency (in the rotating frame of the star) and rotating opposite to the 
stellar spin. This component can excite inertial waves in the convective envelope of the 
star, and the dissipation of inertial waves then leads to a spin-orbit alignment torque 
and a spin-down torque, but not orbital decay. By contrast, for aligned systems, such 
inertial wave excitation is forbidden since the tidal forcing frequency is much larger 
than the stellar spin frequency. We derive a general effective tidal evolution theory for 
misaligned binaries, taking account of different tidal responses and dissipation rates 
for different tidal forcing components. 

Key words: planetary systems - stars: rotation - binaries: close - hydrodynamics - 
waves 



1 INTRODUCTION 

1.1 Tidal Dissipation in Planet-Hosting Stars: A 
Conundrum 

Recent measurements of stellar obliquity in transiting exo- 
planetary systems using the Rossiter-McLaughlin effect have 
shown that a significant fraction of hot Jupiter systems have 
misaligned stellar spin with respect to the planetary an- 
gular momentum axis (e.g., Hebrard et al. 2008; Winn et 
al. 2009,2010,2011; Johnson et al. 2009; Narita et al. 2009; 
Triaud et al. 2010; Pont et al. 2010; Simpson et al. 2011; 
Moutou et al. 2011). This suggests that a large popula- 
tion hot Jupiters are formed through dynamical planet- 
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planet scatterings (e.g., Rasio & Ford 1996; Weidenschilling 
& Marzari 1996; Zhou et al. 2007; Chatterjee et al. 2008; 
Juric & Tremaine 2008) and more importantly, through sec- 
ular interactions between multiple planets or Kozai effect 
induced by a distant companion (e.g., Wu & Murray 2003; 
Fabrycky & Tremaine 2007; Nagasawa et al. 2008; Naoz et 
al. 2011a,b; Wu & Lithwick 2011; Katz et al. 2011), although 
other effects involving star-disc interactions (Lai et al. 2011; 
Foucart & Lai 2011) and the assembly of protoplanetary 
discs (Bate et al. 2010; Thies et al. 2011) may also play a 
role in producing spin-orbit misalignment. 

Recent observations have also revealed an intriguing 
correlation between stellar obliquity O and effective tem- 
perature Tcfi: Winn et al. (2010) found that the misaligned 
systems tend to have hotter parent stars (Tcff ^ 6250 K, cor- 
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responding to stellar mass M^. ^ 1.3Mq), while the systems 
containing cooler stars have small obliquities. This trend 
was also found from a recent analysis of the stellar rotation 
velocities (Schlaufman 2010). Although this Q-Ten correla- 
tion may indicate different planet formation mechanisms for 
stars of different masses, Winn et al. (2010) pointed out 
that Teff = 6250 K corresponds to the temperature be- 
low which stars contain a large sub-surface convective en- 
velope, and suggested that tidal dissipation in these stars 
lead to the damping of their obliquities. Indeed, the facts 
that effective temperature is more closely related to obliq- 
uity than stellar mass and that a few low-mass and long- 
period planets are exceptions to the O — Toff correlation 
(Winn et al. 2010,2011), support the idea of tidal damping 
of spin-orbit misalignment in solar-type stars. Most recently, 
Triaud (2011) found a correlation between the obliquity and 
stellar age, suggesting a tidal alignment timescale of about 
2.5 Gyr. 

It has long been recognized that tidal dissipation in 
planet-hosting stars could play an important role the evolu- 
tion and survival of hot Jupiter systems (e.g., Lin et al. 1996; 
Rasio et al. 1996; Marcy et al. 1997; Sasselov 2003; Dobbs- 
Dixon et al. 2004; Barker & Ogilvie 2009; Jackson et al. 2009; 
Levrard et al. 2009; Leconte et al. 2010; Matsumura et 
al. 2010). The strength of tidal dissipation is usually pa- 
rameterized by a dimensionless quality factor Q* . When the 
orbital mean motion fl is larger than the stellar spin fre- 
quency Cls, the orbital decay timescale (for circular orbits) 
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where Qi = 3Q*/(2fc2) is the reduced tidal quality fac- 
tor, k2 is the Love number, and is the mean density of 
the star and pQ is the solar value. The inspiral time for a 
planet into its host star is (2/13)ta. The observations of hot 
Jupiters with the shortest orbital periods (such as WASP- 
18b, 0.94 d; WASP-19b, 0.79 d; and WASP-43b, 0.81 d; Hel- 
lier et al. 2009,2011; Hebb et al. 2010) suggest Q'^ ^ 10* (see 
Brown et al. 2011). Such a large Q'^, value is consistent with 
recent theoretical works (see Sec. 1.2) on the physics of tidal 
dissipation in planet-hosting solar-type stars (Ogilvie & Lin 
2007; Barker & Ogilvie 2010,2011; Penev & Sasselov 2011). 
On the other hand, in the often-used tidal evolution equa- 
tions (see Sec. 2.5), the damping time for stellar obliquity O 
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where S and L are the stellar spin and orbital angular 
momenta, respectively, Ps is the spin period and the mo- 
ment of inertia of star is kMB?. Thus, for typical param- 
eters of hot Jupiter systems, ~ ta, and a tidal qual- 
ity factor ^ 10* would not cause significant damping 
of O. To put it in another way, a reduction in the stellar 
obliquity is accompanied by a similar amount of orbital de- 
cay, Aa/a ~ (2S/L)AO/ sinO. This poses a severe problem 



to the (otherwise appealing) tidal damping interpretation 
of the stellar obliquity - effective temperature correlation 
(Winn et al. 2010). One way out of this problem is to as- 
sume that the star's convective envelope is weakly coupled to 
its radiative core, thus reducing the obliquity damping time 
for the stellar envelope (Winn et al. 2010). This assumption, 
however, is difficult to substantiate, as fluid instabilities may 
develop in the presence of large differential rotation (espe- 
cially when the directions of rotation vary across the star) 
to quickly couple the rotations of the core and the envelope. 

To recapitulate, there is a conundrum concerning the 
efficiency of tidal dissipation in planet-hosting solar-type 
stars: On the one hand, the survival of hot Jupiters with 
shortest orbital periods and recent theoretical works both 
indicate that stellar tidal dissipation induces only modest 
or negligible orbital decay. On the other hand, the observed 
stellar obliquity - effective temperature correlation suggests 
that tidal dissipation is important in damping stellar obliq- 
uity 

In this paper, we show that tidal damping of spin-orbit 
misalignment can be much more efficient than tidal damp- 
ing of the orbit. In another word, the effective tidal quality 
factor for the former process can be much smaller than the 
latter. This provides a natural resolution to the conundrum 
discussed above. 

1.2 Basic Idea 

Many previous works on tidal evolution in hot Jupiter sys- 
tems (e.g., Rasio et al. 1996; Sasselov 2003; Dobbs-Dixon et 
al. 2004; Barker & Ogilvie 2009; Jackson et al. 2009; Levrard 
et al. 2009; Hansen 2010; Matsumura et al. 2010) were based 
on the weak friction theory of equilibrium tides. This the- 
ory considers large-scale quadrupole distortion of the star, 
and parameterizes tidal dissipation by a dimensionless qual- 
ity factor Qi, or more generally, by a constant tidal lag time 
AtL- The theory was first formulated by Darwin (1880), and 
extensively applied to solar-system bodies (e.g., Goldreich & 
Solter 1966) and stellar binaries (see Zahn 2008 for a review). 
These applications have proved very useful since they pro- 
vide empirical estimates or constraints on the values of 
for various systems. The most general (arbitrary orbital ec- 
centricity and spin-orbit inclination angle) and correct equa- 
tions for tidal evolution based on this theory were derived by 
Alexander (1973), and were also elaborated by others (e.g.. 
Hut 1981; Eggleton et al. 1998; Correia & Laskar 2010). 

Although it is well recognized that the equilibrium tide 
theory is a parameterized theory, with all the physics of tidal 
dissipation hidden in a single parameter Q* or AtL (the Love 
number of the body can be absorbed into the definition of 
these parameters), it is not widely appreciated that the ef- 
fective tidal Q* for different processes (e.g., spin-orbit align- 
ment and orbital decay) can be different. In another word, 
the widely-used tidal evolution equations based on equilib- 
rium tide theory can be incorrect even at the parameterized 
level. Indeed, the conundrum discussed in Sect. 1.1 arises 
because Eq. ^ assumes that the tidal Q* for stellar obliq- 
uity damping is similar to the for orbital decay. In fact, 
this is incorrect, as we explain below. 

There are three channels of tidal dissipations in solar- 
type stars: 

(i) Equilibrium tides. The large-scale quasi-static tidal 
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bulge can be damped by turbulent viscosity in the star's con- 
vective envelope (Zahn 1977,1989). The major uncertainty 
involves how the effective viscosity derived from the mixing- 
length theory, ~ Vth/S (where vt and h are the velocity 
and size of convective eddies, respectively), is reduced when 
the tidal forcing period Ptide is shorter than the convective 
turnover time Tt = h/vt (see Goodman & Oh 1997). Re- 
cent simulations (Penev et al. 2009,2011) suggests that the 
reduction factor is about Ptidc/(27rrt) (for a limited range 
of Ptide), and the corresponding tidal Q^, well exceeds 10^ 
(Penev & Sasselov 2011). An even larger Q* will result if 
the reduction factor (Ptidc/27rrt)^ is used (see Ogilvie & Lin 
2007). 

(ii) Excitation and damping of internal gravity waves 
(Goodman & Dickson 1998; Ogilvie & Lin 2007; Barker & 
Ogilvie 2010,2011). These waves (also called Hough waves 
when modified by rotation) are launched at the bottom of 
the star's convective envelope and propagate toward the stel- 
lar center. If they attain sufficient amplitudes at the center, 
wave breaking will occur; this will produce significant tidal 
dissipation, corresponding to Qi, ~ a few x 10^(P/1 day)*'''^ 
[assuming the orbital period P is much shorter than the 
spin period; see Barker & Ogilvie (2010)]. If the waves 
are reflected coherently at the stellar center (e.g., by a 
small convective core) before nonlinear breaking (Terquem 
et al. 1998), only weak dissipation will result (Q* ^ 10*). 
The latest calculations by Barker & Ogilvie (2010,2011) sug- 
gest that while the nonlinear wave breaking certainly occurs 
for binary stars, it is probably unimportant for exoplanetary 
systems - this would explain the survival of short-period hot 
Jupiters against orbital decay (see Weinberg et al. 2011). 

(iii) Excitation and damping of inertial waves. Recent 
theoretical works on dynamical tides in rotating planets 
(Ogilvie & Lin 2004,2007; Ogilvie 2005,2009; Wu 2005a,b; 
Papaloizou & Ivanov 2005; Goodman & Lackner 2009) and 
stars (Savonije & Papaloizou 1997; Papaloizou & Savoniji 
1997; Savonije & Witte 2002); Ogilvie & Lin 2007) have em- 
phasized the importance of inertial waves driven by Corio- 
lis force. When the tidal forcing frequency (in the rotating 
frame) lj is less than twice the spin frequency {0.s), short- 
wavelength inertial waves can be excited. In particular, when 
these waves are confined to a spherical shell (as in the con- 
vection zone outside the solid core of a giant planet or in 
the convective envelope of a solar- type star), tidal distur- 
bances are concentrated in very narrow regions (called "wave 
attractors") where dissipation takes place (Ogilvie & Lin 
2004; Ogilvie 2009; Goodman & Lackner 2009). It appears 
that this mechanism can explain the tidal Q (~ lO'') of gi- 
ant planets and, when combined with internal gravity wave 
damping [see (ii) above] , can also explain the dissipation re- 
quired for the circularization of stellar binaries. However, 
for solar-type stars in hot Jupiter systems, the inertial wave 
dissipation mechanism is not expected to operate, since the 
tidal frequency iD — 2(il — fis) (assuming circular orbit and 
aligned stellar spin) is larger than 2Qs for typical parameters 
(e.g., P, ~ 10 d and P~ 1 d). 

The main point of our paper concerns inertial wave dis- 
sipation in the parent stars of hot Jupiter systems when 
the stellar spin S is misaligned with the orbital angular 
momentum L. For a circular binary, in the inertial coordi- 
nate system with the Z-axis along L, the tidal potential has 
two components (to the quadrupole order), with frequencies 



LO^i = m'fl (where m' = 0, 2). As seen in the rotating frame 
of the star, the tidal frequencies become 

<imm' = m'Q — mQs, (3) 

with m — 0,±1,±2. For an aligned system, only the m = 
m' = 2 component of the tidal potential is nonzero and in- 
volved in tidal dissipation. (The m = m' = component 
is also nonzero, but it does not transfer energy and angu- 
lar momentum since it is completely static.) For misaligned 
systems, however, all seven tidal components [all combina- 
tions of (m, m') except m = m' = 0; note that for m' — 0, 
the m = — 1(— 2) component is physically identical to the 
m = 1(2) component] contribute to the transfers of tidal 
energy and/or angular momentum. In general, each of the 
7 distinct components will generate tidal disturbance with 
its own quality factor Qmm'- Of particular interest is the 
(m,m') = (±1,0) components. They have ilimm' ~ 
and the angular pattern frequency is a)/m = — f2s (the neg- 
ative sign means it is retrograde with respect to the spin). 
They can generate inertial waves in the convection zone, and 
therefore making a significant or dominant contribution to 
the tidal alignment of spin-orbit inclination. These compo- 
nents, however, do not contribute to orbital decay since they 
are static in the inertial frame and do not transfer energy 
(see Sect. 3). 

In summary, while for aligned hot Jupiter systems, 
the effective stellar tidal Q* governing orbital decay (and 
contributing to orbital circularizatiorij may be quite large 
(^ 10*) because channel (i) and channel (ii) (i.e., equilib- 
rium tides and gravity waves) are ineffective and channel 
(iii) (inertial waves) is forbidden, for misaligned systems, in- 
ertial wave excitation becomes possible, which may provide 
an efficient damping mechanism for the spin-orbit misalign- 
ment. 

The remainder of this paper is organized as follows. In 
Sect. 2 we develop a general effective theory for tidal evo- 
lution of misaligned binaries. It is "general", because the 
theory takes account of the different responses (including 
both equilibrium and dynamical tides) of the star to differ- 
ent frequency components of the tidal potential. It is "ef- 
fective" , because the responses and dissipations of different 
components are treated in a parameterized way. In Sect. 3 
we present the tidal evolution equations due to inertial wave 
dissipation for misaligned hot Jupiter systems. We conclude 
in Sect. 4. 



2 TIDAL EVOLUTION OF MISALIGNED 
BINARIES: A GENERAL EFFECTIVE 
THEORY 

2.1 Tidal Potential 

We consider a star of mass M, radius R and spin Qs (along 
the direction S) orbited by a companion (planet) of mass 
Af'. The orbital semi-major axis is a and the orbital angu- 
lar frequency is il. We allow for general spin-orbit inclina- 
tion angle Q (the angle between the spin angular momentum 
S — SS and the orbital angular momentum L = -LL), but 

^ Tidal dissipation in the planet also contributes to circulariza- 
tion. 
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consider circular orbit for simplicity. In the spherical coor- 
dinate system centered on M with the Z-axis along L, the 
tidal potential produced by M' can be expanded in terms of 
spherical harmonics: 

U{r, t) = -GM' J2 ^^^^ e-^-'^Va^, {OlAl), (4) 



where m' = 0, ±2, with W20 = -(tt/S)^''^ and W2±2 = 
(Stt/IO)^^^. To study the dynamical response of stellar fluid 
to the tidal forcing, we need to express U{r,t) in terms of 
Yim{9, (f>), the spherical harmonic function defined in the in- 
ertial frame centered on M with the z-axis along S. This is 
achieved by the relation 



(5) 



where ©^^/(S) is the Wigner D-matrix of i = 2 (e.g., 
Wybourne 1974), and we have chosen the y-axis along the 
direction S x L. The relevant T>mm' 's are 



2?2±2 = i(l±cose)^ 
V2±i = sine(l ± cose), 

^ ^ . 2^ 
V20 = —r sm 9, 
4 

Dio — sine cose, 

1 2 



and 



m. — m' 



The tidal potential then becomes 

U{v,t) = ~Y, Urnm' rV2„(e, cf,) c' 

mm' 

where 



u„ 



GM' 



-Ur, 
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2.2 Ansatz for Tidal Response 

Each (mm')-component of the tidal potential drives fluid 
perturbation inside the star, which can be specifled by the 
Lagrangian displacement (r, t) and Eulerian density 

perturbation Sp^^i (r, t). In the absence of dissipation, these 
perturbations are proportional to ([/m^//tJo)e""'^~""'"* 
(where = ^GM/R-^ is the dynamical frequency of the 
star), exactly in phase with the tidal potential. When fluid 
dissipation is present, there will be phase shift between the 
fluid perturbation and the tidal potential. This phase shift, 
in general, depends on m, the forcing frequency (in the rotat- 
ing frame) oJmm' , as well as the intrinsic property (including 
the rotation rate) of the star. We write this phase shift as 

^mm' ^mm'^mm' ^ (1^) 

where t^^i is the "lag time" associated with the {mm')- 
perturbation. Thus, without loss of generality, we can write 
(mm')-component of the fluid displacement as 



-imm' (r) exp(-im'm -I- iA„ 



(16) 



— 5pmm'(r) exp(-im'J7f + iA^^;), (17) 



The corresponding density perturbation is 

'5Pmm'(r,t) = — 



with 

5pw =-V-(p|„„,)- (18) 

Note that ^^^/(r) and (5pmm'(r) are proportional to e™*^; 
except for this factor, they are real functions. Also note that 



(19) 



with = W2m''Dmm'{Q)- Thus, (r) and 5p™„/(r) 

specify the amplitudes (and shapes) the fluid perturbations 
after factoring out the "equilibrium" dimensionless tidal dis- 
tortion, {M' /M){R/af. 

Note that although our ansatz of the tidal responses is 
motivated by the weak friction theory of equilibrium tides, 
it is actually quite general, provided that one keeps in mind 
that the lag time t„im' depends on m, uj^m' and the intrinsic 
property of the star. 



Note that when expressed in terms of (j>r (the azimuthal an- 
gle in the rotating frame of the star), each term in Eq. H12| 



has the dependence 6*""^'-+™"=*-™''^*. Thus, the tidal po- 
tential is composed of various (mm' )-components, each with 
forcing frequency m'Q. in the inertial frame. In the frame 
corotating with the star, the forcing frequency of the (mm')- 
component is Umm' ~ m'Q, — mQ.s [Eq. ((3])], with the pattern 
rotation frequency 



— 1 J7 - n« 
m 



(14) 



(Obviously, the pattern frequency has no meaning for 
the axisymmetric m = components.) Note that 
physically, there are 7 distinct components: (m, m!) = 
(0, 2), (±1, 2), (±2, 2), (1, 0), (2, 0) (see also Barker & Ogilvie 
2009). The (0, 0) component is static and does not contribute 
to tidal dissipation. The (— m, — m')-component is physically 
identical to the (m, m')-component. 



2.3 Tidal Torque and Energy Transfer Rate 

The tidal torque on the star is 

ifx <5/9(r , t) r X [- V(7* (r , t)] , 



(20) 



and the tidal energy transfer rate (from the orbit to the star) 
is 

E = jd?x p(r) . h Vt/- (r , t)] . (21) 

Using Eqs. (|12p and (|17p . the z-component (along the spin 
axis) of the tidal torque reduces to 



where 



(22) 



(23) 
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and Krnm' IS the "Love coefficient" 



^ - /d\5pw(r)rV;„(e,<^). (24) 



Similaly, the energy transfer rate is given by 

E = TqD. ^ Ul^rn'm! timm' sin A„^/ . (25) 

mm' 

Note that the (mm' )-component of energy transfer and tidal 
torque satisfies 



(-^)mm' (Tz^mm' 

m 



(26) 



This is expected since the pattern rotation rate (in the in- 
ertial frame) of the tidal force is m'Q,/m. 

Equations (|22p and H25p can be worked out explicitly. 
Assuming jAm^/l <JC 1, we have 



Tz = 



^^To [(1 + ce)\n - ^s)t22 

-se(l-ce)''(2n + n,)r_i2 
-(l-C9)^(n + n,)r_22] 

— 3(H^2o)^rO (Sg T20 + sl Ce no) 



(27) 



and 



£ = -^^^Sy^Tof^ [(1 + ce)*(n - f7,)r22 

+2s^(l + ce f{2rL - VLs)ti2 + 6s^ J7ro2 
+2s^(l-ce)'(2!^ + f7,)T_i2 

+ (l-C9)^(n + !^,)r_22], (28) 
where we have defined 

and cg = COS0, sg = sinG, and we have used the identity 

The perpendicular component of the tidal torque T 
is less straightforward to evaluate. The y-component does 
not depend on tidal dissipation (for |A^^/| <C 1) and con- 
tributes to the spin precession. Since Ty does change the 
spin-orbit elements, we will not consider it further. The x- 
component of T is 



rmnm 5 



(30) 



where 



tan t 



-Yoff, cos( 



>)■ 



(31) 



Since 5p^^i{v) oc e™**, Krnfnm' is nonzero only for m 
m ± 1. If we define 



(32) 



then Eq. (|30p can be written as 

jnjn' 

(33) 



A direct calculation shows that the coefficients k 

n 

related to the Love coefficients [Eq. (|24p ] by 



^2m' ^ ^2m' , 



— K — 2n 



^3/2 Kl^ 



/372ko„ 



+ 



"--im' ^ V 3/2 «;_!„/, — 
Thus Ta; can be worked out explicitly: 



(34) 



Tx = 



(W2: 



\2 r 

^ro[se(l + ce)^(!^-n,)r22 



+se(l + ce)^(2 - C9)(2!:! - n,)ri2 

+ 6 Sg Q, T02 

+sg{l - cgf{2 + ce)(2n + n«)r_i2 
+se(l-ce)-^(!:! + n,)r_22] 
+3(^20)^ To (s^ Co r20 + Sg 4 no) 



(35) 



2.4 Tidal Evolution Equations 

Given the tidal torque Tz and energy transfer rate E, the 
tidal evolution equations for the stellar spin Q,s and the or- 
bital semi-major axis a are 



2aE 
' GMM' 



(36) 



where I is the moment of the inertia of the star. The spin- 
orbit misalignment angle Q is given by cos = S ■ 1,/{SL), 
where S = IQs and L — fj,a^fl. Using S — — L — T, we find 

A Nx Nx Nz 

9 = — -; — COS0 H — — sm9. (37) 

J -L L/ 

Note that the rate of change for the magnitude of the orbital 
angular momentum satisfies 



L = -Tz cose - Tx sine = -| = 



(38) 



This can be checked directly using Eqs. (p7|) . (|28|) and (|35|) . 
Thus a circular orbit will remain circular, as it should be. 

These evolution equations for Cls, a and Q [with 
Tz, E, Tx given by Eqs. §5^, respectively] 

are the most general tidal equations for circular binaries. 
They depend on 7 independent "reduced" tidal lag times 
Tmm' = i-mm' i^mm' , Corresponding to the 7 independent tidal 
forcing components. In general, each Tmm' depends on m, 
ujrnm' ^iid the iutriusic property (including Q,s) of the star. 



2.5 Special Case: Weak Friction Theory of 
Equilibrium Tide 

When Tmm' = T are the same for all seven tidal forcing 
components, we find 

E ='^TQQ.(n-Q.sCosQ^r, (39) 

Tz = ^To [2f7 cose - (1 + cos^e) Os] r, (40) 

Tx = ^Tq sine {20. - Q.S cose^ r. (41) 
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These give 



a 
a 

e : 
where 



1-^cose 



L 



ta \2S 

1 f L_ 

2S 



cos O 



sin O 



|)(i + cos^e) 



cosG 



and we have defined 



3fc2 

Q 

4% 



k2AtL = -^T, Q= (2nAtL) \ 

5 ^ ' 



(42) 
(43) 
(44) 

(45) 
(46) 



Here fc2 and /SIl have the usual meanings as in the equili- 
birum tide theory: k2 is the Love number and AtL is the 
tidal lag time [which is related to the viscous time tvis 
by Ati, = l/(a;otvis)-] These equations agree with those 
given in Alexander (1973) and others (Hut 1980; Eggleton 
et al. 1998). 



3 TIDAL TORQUE DUE TO INERTIAL WAVE 
DISSIPATION 

The convective envelope of a solar-type star supports iner- 
tial waves, which are driven entirely by Coriolis force. The 
frequency uj (in the rotating frame) of an inertial wave is 
related to its (local) wavenumber vector k by the dispersion 
relation (e.g., Greenspan 1968) 



(2n, ■ k/|kj)' 



(47) 



Thus inertial waves exist only when \uj\ < 2Qs. 

In hot Jupiter systems with f2 S> Os, the only tidal 
forcing component that is capable of exciting inertial waves 
in the star is (m,m') = (1,0), with the forcing frequency 
(in the rotating frame) uj = —Qa- [The (— 1, 0)-component 
has UJ = Qa and is physically identical.] The tidal torque 
associated with this component can be read off directly from 
Eqs. (pTj) and p5|l : 



(48) 
(49) 



iT^)io = -^ToQaTw (sine cose) ^ , 

(r,)io = ^To^aTio sinecos^e, 



with To = G{M'/a^fR^. Similar to Eq. ((46|, we define the 
relevant tidal Love number fcio, lag time Atio and quality 
factor Qio: 



47T" 

fcioAtio = — Tio, Qio = (risAtio)" 
5 

Using Eqs. (O and (jSTll, we find 



^) =_^(sinecose)^ 

"s / iQ IslO 

''e) = — 3_ sine cos^e (cose + ^ 

V / 10 talO V L 



(50) 

(51) 
(52) 



where 

1 _ 37rrioTo 

talO 5/ 



(53) 



40io V M , 

For Al = Af* and Al' — Alp, the corresponding timescale is 
io/fcio\ f AL 



t..o = 4.3(- 
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\wnipj \pQ 



(54) 



[cf. Eq. ©]. 

Since the (1, 0)-component of the tidal potential is static 
in the inertial frame, the associated energy transfer rate E 
associated with this component is zero, giving 



0. 



(55) 



This does not imply zero tidal dissipation. In fact, the tidal 
energy dissipation rate equals the energy transfer rate in the 
rotating frame, 



Er = -QaT^ > 0. 



(56) 



This is exactly balanced by {d/dt){S^ /2I) = ^aS = fi^Tz < 
0, so that E = Er + ^aT^ = 0. Also note that the (1, 0)-tidal 
force does not change the magnitude of the orbital angular 
momentum, L = —T^ cose — Ta; sine = 0. So a circular orbit 
will remain circular. 

Two interesting features of Eq. (|52[) are worth noting: 
(i) e = when e = 90°. Thus, if the (l,0)-tidal component 
is the only tidal force operating in the system, there could 
be many systems with the spin-orbit misalignment angle 
stalled around 90°. (ii) 6 > when cose < -S/L. This 
implies that a retrograde system (e > vr/2) may evolve to- 
ward anti-alignment. Of course, other tidal components can 
also contribute to 6 (Sect .2) and will weaken these features. 
But with enough statistics of misaligned systems, it may be 
possible to test these or constrain fcio/Qio. 

A detailed calculation of Qio/fcio is beyond the scope 
of this paper. Previous works of tidal dissipation in rotat- 
ing solar-type stars have focused on the m = 2 tide for 
aligned binaries (Savonije & Papaloizou 1997; Papaloizou & 
Savoniji 1997; Savonije & Witte 2002; Ogilvie & Lin 2007) 
and demonstrated the importance of inertial waves. In the 
numerical study by Ogilvie & Lin (2007), inertial waves af- 
fect tidal dissipation in two ways: (i) short-wavelength in- 
ertial waves in the convection zone are damped by turbu- 
lent viscosity; (ii) Inertial waves influences the excitation of 
gravity (Hough) waves in the radiative zone. Regarding (ii), 
they assumed that the inward propagating Hough waves are 
damped near the stellar center; this is appropriate for binary 
stars, but may lead to an over-estimate of stellar tidal dissi- 
pation in hot Jupiter systems (Barker & Ogilvie 2010,2011; 
see Sect. 1.2). Ogilvie & Lin (2007) found that when iner- 
tial waves are excited {\uj\ < 2fla), the tidal dissipation rate 
depends on uj in an erratic manner. On average, the energy 
dissipation rate in the convection zone is significantly in- 
creased (by 1-3 orders of magntitude, depending on the spin 
period) compared to equilibrium tides. For example. Figure 
3-6 of Ogilvie & Lin (2007) show that for a solar-type star 
a,t UJ — —Qa, the tidal quality factor associated with iner- 
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tial wave dissipation is ~ 6 x 10^ for Ps — 10 days and 
Q'^c^5x for P, = 3 days. 

Barker & Ogilvie (2009) reported the result of a calcula- 
tion of the I = 2, m — 1 tidal dissipation for a specific F-type 
stellar model (appropriate for the hot Jupiter system XO-3) 
using the numerical method of Ogilvie & Lin (2007). In- 
triguingly, they found that at iIi/Qs = —1, tidal dissipation 
is significantly enhanced, with the effective tidal Q ~ 10® 
(see their Fig. 7). They suggested that this prominent fea- 
ture arises from resonant excitation of the I = m — 1 Rossby 
mode. More systematic studies on the {m,m') — (1,0) tide 
for a range of stellar models (with different sizes of the con- 
vective envelope) would be useful. 



Q parameters may be adopted. They are 



4 CONCLUSION 

The main point of this paper is that in close-in exoplane- 
tary systems, when the stellar spin axis is misaligned with 
the orbital angular momentum axis, a new tidal dissipation 
channel opens up. This channel involves the excitation of in- 
ertial waves in the stellar convection zone, and is forbidden 
for aligned systems. Thus, tidal damping of spin-orbit mis- 
alignment can be more efficient than orbital decay. This may 
explain the stellar obliquity - effective temperature correla- 
tion observed by Winn et al. (2010) and the obliquity - age 
correlation noted by Triaud (2011), while still being consis- 
tent with the survival of hot Jupiters with very short orbital 
period (see Sect. 1.1). 

On a more general level, this paper highlights the im- 
portance of treating tidal dissipation as being dependent 
on the tidal forcing frequency, the strength of tidal poten- 
tial, and the tidal processes involved (e.g., orbital decay vs. 
spin-orbit alignment), in contrast to the equilibrium tide 
equations widely used in many applications and empirical 
works (see references in Sect. 1.2). Indeed, various studies of 
the physical mechanisms of tidal dissipation (e.g., Ogilvie & 
Lin 2004,2007; Goodman & Lackner 2009; Barker & Ogilvie 
2010,2011; see also Zahn 2008 for a review of earlier works) 
have already made this point clear. While it is recognized 
that the weak friction theory of equilibrium tides (Darwin 
1880; Goldreich & Soter 1966; Alexander 1973; Hut 1981; 
Eggleton et al. 1998) is a parameterized theory (with the 
tidal Q or lag time being the single parameter), we have 
shown in this paper that even at the parameterized level, 
its equations are sometimes inadequate or misleading, since 
different tidal processes (e.g., orbital decay vs. spin-orbit 
alignment) may involve very different tidal dissipation mech- 
anisms. 

Given the complicated nature of tidal dissipation, a pa- 
rameterized "effective" theory of tidal evolution remains use- 
ful. We have derived such an effective theory in this paper 
(see Sect. 2). For misaligned circular binaries, there are seven 
independent tidal quality factors Qmm' or lag times Tmm' in 
the theory, related to the tidal responses for different Fourier 
components of the tidal potential. Obviously, such a system 
of equations with 7 parameters (Sect. 2.4) is not convenient 
to use in real applications. For hot Jupiter systems, we sug- 
gest that a reduced system of equations involving two tidal 



a 
a 

sis 



+ 



cq \ / 10 

e = (e)eq + (e)io - (e)io,oq 



10, cq 



(57) 

(58) 
(59) 



where (a/a)cq, {Cls/fls)cq and (0)cq are given by Eqs. (|42p - 
(|44|l and characterized by the parameter k^jQ, while 
(f2s/f2s)io and (O)io are given by Eqs. (|51|l - (|52| l and char- 
acterized by the parameter fcio/Qio. The expressions for 
(f2s/f2s)io,oq and (O)io,oq are the same as Eqs. (f5T]) - ([52)) . 
except that the the parameter fcio/Qio should be replaced 
by fe/Q. These equations (extended to eccentric orbits), 
combined with similar equations for planetary tides, can be 
used to assess and constrain the effects of tidal evolution in 
close-in exoplanetary systems. 
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